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We show that classical density functional theory can be based on the constrained search method 
[M. Levy, Proc. Natl. Acad. Sci. 76, 6062 (1979)]. From the Gibbs inequality one first derives a 
variational principle for the grand potential as a functional of a trial many-body distribution. This 
functional is minimized in two stages. The first step consists of a constrained search of all many- 
body distributions that generate a given one-body density. The result can be split into internal 
and external contributions to the total grand potential. In contrast to the original approach by 
Mermin and Evans, here the intrinsic Helmholtz free energy functional is defined by an explicit 
expression that does not refer to an external potential in order to generate the given one-body 
density. The second step consists of minimizing with respect to the one-body density. We show that 
this framework can be applied in a straightforward way to the canonical ensemble. 
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I. INTRODUCTION 

The variational principle of density functional theory 
(DFT) was originally formulated for ground-state prop- 
erties of quantum systems by Hohenberg and Kohn in 
1964 The extension to non-zero temperatures was 
performed by Mermin in the following year Q, here still 
formulated for quantum systems. The application to the 
statistical mechanics of classical systems, i.e. the devel- 
opment of classical DFT, was initiated about a decade 
later through the work of Ebner, Saam, and Stroud 

a- 

5] . The generality of the framework was fully realized by 
Evans [6j]. His 1979 article continues to be the standard 
reference on the subject; there are more recent review 
and textbook [l(| presentations. 

The Hohenberg-Kohn theorem applies to one-particle 
density distributions that correspond to a particular ex- 
ternal (one-body) potential v in the Hamiltonian, in 
which the kinetic and internal interaction are those of the 
true system [l], Q ■ One refers to v-representability of the 
one-particle density, i.e. the condition that a one-particle 
density is generated by some external potential v. How- 
ever, it was realized, already in the original Hohenberg- 
Kohn paper, that w-representability is n ot g uaranteed for 
an arbitrarily chosen density p [l], LLl|, LL2[ • One argues 
that this does not pose any problems in the practical ap- 
plications of DFT to quantum systems [13[ . In the devel- 
opment of the theory, it turned out that there indeed ex- 
ist non-u-representable densities, i.e. one-body densities 
that are not associated to any ground-state wave func- 
tion [3]. The original Hohenberg-Kohn theorem does 
not apply to these. 

In 1979 Levy introduced an alternative foundation of 
DFT for quantum systems, based on a constrained, two- 



stage search [ll| . Here a weaker condition, known as N- 
representability, is used, where the density distribution 
may be directly obtained from some anti-symmetric TV- 
body wave function, although an external potential that 
generates this wave function need not exist [HI EH EH ■ 
One defines an exchange-correlation functional that de- 
mands searching all wave functions that return the fixed 
(trial) one-body density. The latter need not be v- 
representable. Subsequently, a method similar to Levy's 
was proposed by Lieb called the generalized Legen- 
dre transform [1 21 ] . Instead of searching all wave func- 
tions, the functional searches all possible external poten- 
tials that correspond to a fixed density. Kohn adopted 
the constrained search for his Nobel lecture [Hj], and it 
is viewed as an important theoretical contribution to the 
foundation of DFT for electronic structure. Practical ap- 
plications of constrained search functionals are of ongo- 
ing research interest, see e.g. [IH, E3- Levy gives a brief 
historic account of the development of his ideas in Ref. 

Given the significance of Levy's and Lieb's methods 
for electronic structure, it is somewhat surprising that 
there are very few studies that point to the use of these in 
classical systems. One example is the work by Weeks [19j , 
where the w-representability of the one-body density in 
some finite region of space is investigated through the 
Gibbs inequality. Although Weeks cites Levy's original 
paper (Tlj . and makes a remark that his formulation is 
in spirit similar to that of Levy, it seems that his method 
is related more to Lieb's generalized Legendre transform 
method. Earlier work has been carried out in order to 
investigate the existence of an external potential that is 
associated to a given equilibrium one-body density pol |. 
It is concluded that there is such an external potential 
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that produces any given density for any (classical) system 
without hard core interaction. Although one might guess 
from general arguments that the constrained search can 
be applied to classical DFT, to the best of our knowledge, 
this procedure has not been spelled out explicitly in the 
literature. 



In the present article we show how to formulate the 
variational principle of classical DFT based on Levy's 
constrained search method. This alternative can provide 
further insights into the foundation of classical DFT. In 
particular the intrinsic free energy functional is defined 
here without implicit reference to an external potential 
v. A more relaxed condition for p, similar to that of N- 
representability, is imposed. Here, the one-body density 
is only required to be obtained from an arbitrarily chosen 
many-body probability distribution /. We refer to this 
condition as /-representability of a given p. While dis- 
tinguishing between the different type of representability 
in practical DFT calculations seems unnecessary, we find 
the discrimination very useful for conceptual purposes 
and hence point out throughout the manuscript which 
of the conventions is followed in the reasoning. We also 
show that Levy's method can be applied in a straight- 
forward way to the canonical ensemble. There is consid- 
erable current interest in the theoretical description of 
the behaviour of small systems, where the grand and the 
canonical ensembles are inequivalent in general, and the 
later might model certain (experimental) realizations of 
strongly confined systems more closely. Several relatively 
recent contributions address the problem of formulating 
DFT in the canonical ensemble [2ll - l26| . The authors of 
these papers consider the important problem of how to 
obtain DFT approximations that make computations in 
the canonical ensemble feasible. Our present article has 
a much lower goal: We are only concerned with formu- 
lating the variational principle in an alternative way. 



This article is organized as follows. We start by defin- 
ing the grand potential as a functional of the many-body 
probability distribution in Sec. |TTJ This is a necessary 
step and our presentation follows [l(J. In Sec. Hill 

we give a brief overview of the standard proof of DFT, ex- 
pressing the free energy as a functional of the one-body 
density based on a one-to-one correspondence between 
the one-body density and the external potential. The 
full derivation is widely known and can be found in nu- 
merous references @-S E3, H3l ■ We proceed, in Sec. HV1 
by formulating the intrinsic free energy functional via the 
constrained search method; our presentation is similar to 
Levy's original work Our central result is the def- 

inition (|22[) of the intrinsic free energy functional, with- 
out reference to an external potential. We summarize 
the essence of Levy's argument as a double minimization 
[lH in Sec. |V] In Sec. I VII we apply this to the canonical 
ensemble and we conclude in Sec. IVII1 



II. GRAND POTENTIAL FUNCTIONAL OF 
THE MANY-BODY DISTRIBUTION 

In the grand canonical ensemble of a system of classical 
particles, the equilibrium probability distribution for N 
particles at temperature T is assumed to exist and to be 
given by 

fo = Z- 1 exp(-f3(H N -pN)), (1) 

where Hn is the Hamiltonian of N particles, p is the 
chemical potential, and j3 = l/(fcsT), with ks being the 
Boltzmann constant. The normalization constant is the 
grand canonical partition sum 

S = Tr cl e X p(-/3(ff J v-/ziV)), (2) 

where Tr c i represents the classical trace, i.e. the sum over 
total particle number and integral over all degrees of free- 
dom 

N=0 ' J J 

where h is the Planck constant, ri, . . . ,rjy are the po- 
sition coordinates and pi , . . . , pn are the momenta of 
particles 1, . . . , N. 

One introduces the grand potential as a functional of 
the many-body probability distribution, 

V[f] = Tr cl / (H N - piV + /T 1 In /) , (4) 

where / is a variable trial probability distribution that 
satisfies the normalization condition 

Tr cl /=1. (5) 

Note that / as an argument of functional (|4|) can be quite 
general and need not be linked to an external potential 
at this stage. Inserting the equilibrium probability dis- 
tribution |T]) into ((3]) one obtains 

fi[/o] - Tt cl fo (HN-fiN + 13- 1 hi f ) 

= Tr cl / (H N — piV + /T 1 [-InS - p(H N - piV)]) 

= -/T^nS^o, (6) 

where Qq is the equilibrium grand potential. An im- 
portant property of functional Jl} is that it satisfies the 
variational principle 

fi[/]>n[/o], f*h, (?) 

which may be proven using the Gibbs-Bogoliubov rela- 
tion as follows. First, from JU and ©, fl[f] of Eq. (g| 
can be written as 

Q[f]=n[fo]+r 1 Tr cl fbx(-£j . (8) 
According to the Gibbs inequality @, H3] , 
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and hence 



Tr cl /ln( -f ) >Tr cl (/-/ ) 



(10) 



Since / and fo are normalized, i.e. satisfy ([5]), the RHS 
of the inequality above vanishes, and 



/T^d/lnf^-) >0. 

Jo, 



(11) 



Thus the second term on the RHS of ([8]) is positive, and 
the inequality ([7]) follows. 

For classical particles the Hamiltonian may be re- 
stricted to the form 



A' 



N 



2m 

i=l i=l 

where the first term is the total kinetic energy, with the 
squared momentum pf = p, ■ of the i-th particle, U is 
the interatomic potential between the particles, v is an 
(arbitrary) external one-body potential and m is the par- 
ticle mass. The equilibrium one-body density at position 
r is given as a configurational average 



p (r) = Tr cl / 0/ 5(r), 



(13) 



where the density operator for TV particles is defined as 



(14) 



i=i 



The functional form (j4]) was originally introduced by 
Mermin Q for (finite temperature) quantum systems, 
where the grand potential is a functional of a (trial) den- 
sity matrix. The variational principle ([7]) will be used in 
Sees. IIIII and IIVI below, where we present two alternative 
derivations of the intrinsic free energy as a functional of 
the one-body density. 



III. MERMIN-EVANS DERIVATION OF THE 
FREE ENERGY FUNCTIONAL 

Evans gave a formal proof that the intrinsic free en- 
ergy of a system of classical particles is a functional of 
the one-body density [7J . Here we briefly lay out his chain 
of arguments. The many-body distribution fo as given in 
((T|) is a functional of the external potential v through Eq. 
(fT2"| . and therefore po is a functional of v via (fT3"]) . This, 
in principle, requires solution of the many-body problem 
and the dependence is in accordance with physical intu- 
ition, i.e. it is the action of v that generates the shape of 
the density profiles po- 

However, the more useful result that can be de- 
duced p}, is that once the interatomic interaction po- 
tential U is given, fo is a functional of po- The proof 
of this statement rests on reductio ad absurdum [l(| , 



where for a given interaction potential U , v is uniquely 
determined by po- The resultant v then determines fo 
via ([I]) and ([12")) . Hence, fo is a functional of po. 

An important consequence in this reasoning is that for 
given interaction potential U, 



N 



T\p] = Tr cl / ( Y, P~ + U + /T 1 In /, 



(15) 



is a unique functional of the (trial) one-body density p. 
Here, the dependence of fo on the external potential, v, 
is now only implicit through the one-body density, p. We 
will comment on this sequence of dependences in the con- 
clusions, after having laid out Levy's alternative method 
to define a free energy functional in Sec. IIVI Further- 
more, using a Legendre transform, the grand potential 
functional is obtained for a given external potential as 



n v [p]=J r [p} + J dr(«(r)- M )p(r) 



(16) 



The functional f2 u [p] returns its minimum value if p = po, 
i.e. if the trial density is the true equilibrium one-body 
density of the system under the influence of v. The value 
is the grand potential fto- The existence of the min- 
imum value of Sl v [p\ may be proven by considering an- 
other equilibrium density p' associated with a probability 
distribution /' of unit trace, such that 



= Tr cl / / ( J ff JV -/iiV + /r 1 ln/') 
= T[p']+ J &v{v{v)-p)p>{v) 

= ttv[p'], 



(17) 



where 



JY 



(18) 



However, it is known from Eq. ([7]) that > fi[/o], for 
f ^ fo, thus it follows that 



(19) 



In other words, the correct equilibrium one-body density, 
po, minimizes £l v [p] over all density functions that can be 
associated with a potential v. 

This important result may be stated as a functional 
derivative 



Sfl v [p] 



5p{r) 



= 0, 



and 



fl v [p ] = ^o- 



(20) 



(21) 



To conclude, the formal argument for the definition 
dT51) of the intrinsic free energy functional, F[p\ , is based 
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on u-representability of the one-body density. A v- 
representable p is one which is associated with a proba- 
bility distribution, /, of the given Hamiltonian Hn with 
external potential v (ill . [l3j j . This condition was origi- 
nally introduced for quantum systems, and is implicit in 
the current approach. It is used to prove the chain of de- 
pendency outlined above (|15|) . and confirmed for a large 
class of (classical) systems in |20j . 



IV. FREE ENERGY FUNCTIONAL VIA LEVY'S 
CONSTRAINED SEARCH METHOD 

Here we show how one may alternatively define a free 
energy functional via Levy's method. This is based 
on the weaker condition of /-representability, where 
trial density fields p need not necessarily be associated 
with some external potential. We define the intrinsic 
Hclmholtz free energy functional as 



J" L [p] = min 



JV 



T ^(ES7 +c/ +^ lln / 



2m 



, (22) 



where the minimization searches all probability distribu- 
tions /, that are normalized according to ([5]), and that 
yield the fixed trial one-body density p via 



p(r) = Tr cl /p(r). 



(23) 



The notation / — > p in (|22|) indicates the relationship 
(|2"3")l . Note that i) in general there will be many different 
forms of / that yield the same p, and ii) no further con- 
ditions on / are imposed, apart from its normalization. 
In particular, the form of / need not be of Boltzmann- 
type containing the interaction potential U (as was the 
case in Sec. IIIip . Hence p need only be /-representable, 
but not necessarily u-representable. ^l[p] returns a min- 
imum value by choosing the probability distribution that 
minimizes the term in brackets in (|22|) . Note that the 
functional form of this term is formally equivalent to (|15[) 
and that it is a sum of contributions due to kinetic en- 
ergy, internal interaction energy U, and (negative) en- 
tropy /cs/ln/ multiplied by T. 

The grand potential functional for a given external po- 
tential is then 



n L [p]=.F L [p] + J dv(v(v)-fx)p(r) 



(24) 



This functional possesses two important properties, i) At 
the equilibrium density it yields the equilibrium grand 
potential 



^l[po] = 



where po is given by (|13p and f^o by 
constitutes the minimum such that 

n L [p] > n . 



(25) 



This value also 



In order to prove ([25]) and (|26|) , we introduce additional 
notation. Let / min be the probability distribution that 
satisfies the RHS of Eq. Then it follows that 

/ N 2 \ 
JxM = Tr cl /£ in jX^^L + U + p- 1 In /£.J , (27) 

and for the case of the equilibrium density 

Fl\Po] = Ttd/Sn ( E %t + U + ^ 1 ln/Zl ) ■ (28) 

First we proof the inequality (|26[) . By its very defini- 
tion (f2~4")l. the LHS of ([2"B)) may be rearranged into 



dr (v{r) -/i)p(r) + J" L [p] 
dr (v(r) - p) p(r) 



N 



+Trci/^ n (E£ + C/ + ^ lln /n li n) (29) 



= Tr cl / min (^-piV + r'ln/^). 
But according to the inequality (|7|l. 

Trci/f nin (H N -pN + /T 1 In / min ) > Sh- (30) 



Thus, combining (|29j) and (|30|) . the inequality (|26l) is 
recovered. In order to prove ([25jl . it is obvious from 
that 

Trci/a (Hjv - + r 1 In > n„, (31) 
or, recalling ([6]), 

TW^^-M^V + r'ln/a) 

> Tr c i/ (ffjv - /*JV + /T 1 In / ) . (32) 

fPO 

min 

hence from 



But /^° in and /o generate the same one-body density po, 



dr(u(r) -p) p (r) 

{ N 2 \ 
+Tr cl /^ n f £ + [/ + r l l n j > 

dr (v(r) -p) Po(r) 

/ w 2 \ 
+Trci/o Ei + ^ + ^^o , 



(33) 



we obtain 



(26) 



/ N 2 \ 

TTcl/^ f £ + ^ + Z 3 " 1 ln /min J > 



(34) 
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However, by the very definition of the following 

inequality should also hold: 



JY 



2m 



N 



Tr ^(E£ + c/ +^ lln /« 



(35) 



The above two inequalities hold simultaneously, if and 
only if equality is attained, 



Po 
min 



\i=l 

/ N 2 

= Tr cl / [Y i p- + U + 0- 1 hifo 



(36) 



Inserting ([2"8|1 into (j3"6"]l yields 



AT 



Fl [po] = Tr cl / ( XI f| + + r 1 In /< 



(37) 



Furthermore, as 

fio = Trd/o^AT-M^V + ^Mn/o) 
= / drf>(r) -/i)po(r) 



JV 



T Wo(E|^ + t/ +^ lln /« 



(38) 



inserting (j3"7| into (|3"8"]l returns (f^5|) . which completes the 
proof. Equation (|37j) implies that if p is u-representable, 
then Fl[p] = F[p]. Moreover, /q = /^° in means that 
/o may be obtained directly from po even if v is un- 
known: find the probability distribution which yields po 
and which minimizes (|22[) . 

Finally, the inequality (f26| implies that the functional 
derivative of the grand potential functional vanishes at 
equilibrium, 



5p(r) 



= 0. 



(39) 



For completeness we mention that it is convenient to 
split Fl[p] into two terms, viz. the ideal and excess free 
energy functionals, .Fid I/?] and F exc [p]j respectively, such 
that 



F oxc [p] = F L [p] - F id [p], 



(40) 



where the free energy of the ideal gas (with no interaction 
potential present, U — 0) is given by 



FidW = r 1 J drp(r) (in (A 3 p(r)) - l) 



(41) 



where A = {h 2 (5 / (2rmr)) 1 " . Thermodynamics enters by 
realizing that Fl[po] is the 'intrinsic' Hclmholtz free en- 
ergy of the system, such that the total free energy is the 
sum of internal and external contributions, 



F L [p ] + J drp (i>(r). 



V. TWO-STAGE MINIMIZATION 



(42) 



The essence of the derivation presented in Sec. IIVI is 
a double minimization of the grand potential functional 
((4]) of the many-body distribution. In the following, we 
spell this out more explicitly. From Sec. [II] we know that 

n = mm Tr cl / (H N - pN + /T 1 In /) . (43) 

We decompose the RHS into a double minimization 

fi = min min Tr cl / (H N - pN + /T 1 In /) , (44) 

where the inner minimization is a search under the con- 
traint the / generates p (via relationship (|2"3"10 . For 
Hamiltonians of the form (fT2")) the above can be written 
as 

/ \ p 2 

fin = min min Tr r i f > — - — \- U 
p f^P J 2m 

N \ 

+ 5>(r J )-pJV + /T 1 ln/ . (45) 
»=i / 

In the expression above 

Tr cl / (j2 «( r *) ~nA = J dr (v(r) - p) p(r), (46) 
because / — > p. So we may re-write (|4"5j) as 

tt = min jy dr (v(r) - p) p(r) (47) 



N 



■minTr cl / V + [/ + /T 1 In / 
f^ P \^ 2m 



or 



= mm [J dr (v(r) - p) p(r) + F L [p] 



(48) 



where Fl[p] is given by (|22j) . Clearly (j48|) is equivalent 
to US} and (US]). 



VI. DFT IN THE CANONICAL ENSEMBLE 

One benefit of Levy's method is that it allows straight- 
forward generalization to the canonical ensemble, as we 
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demonstrate in the following. In the canonical ensem- 
ble (i.e. for fixed number of particles, N) the equilibrium 
many-body distribution functions is 

f Nfi = Z - 1 exp(-/3H N ), (49) 

where the canonical partition sum is 

Z Q = Tr N eM-/3H N ), (50) 

with the canonical trace 
1 



Tr 



dri . . . dr N / dpi ... dp 



(51) 



h 3N Nl 

In analogy to (U]) we define the functional 

F[f N ] = TtnIn (H n + /T 1 In f N ) , (52) 

where /jv is an arbitrary TV-body distribution that sat- 
isfies Tin/n = 1. It is easy to show that the (total) 
Helmholtz free energy Fq — — /3 IxiZq is obtained by 
inserting the equilibrium distribution (|49p into the func- 
tional (|52|) . hence 



Fq — F[/n,o\- 



(53) 



Reasoning based on the Gibbs-Bogoliubov inequality, 
completely analogous to the arguments presented in 
Sec. in yields 

F = min F[f N }. (54) 

In 

We decompose this into a double minimization 

F a = min min F[f N ], (55) 

where the canonical one-body density distribution that 
is generated by /at is 



p N {v) = Tr N f N p(r), 



(56) 



with the density operator p(r) defined by (fT4]). Clearly 
the density defined in this way satisfies J drp]y(r) = N, 
and there are no fluctuations in the total number of parti- 
cles. For Hamiltonians of the form (fT2j) . Eq. (|55|) becomes 



N 9 

Pi 



F = min min Tr N fN V" „ 

pn f N ~*PN \ f— f 2m 



N \ 

+ C/ + ^«(r ! )+/3- 1 ln/AT 
i=i / 



(57) 



In the above expression 

N 



TrAr/Ar ^w(ri) = J drv(r)p N (r) 



(58) 



because /at — > pat via (|56j) . Hence 



mm 

PN 



drv(r)p N (r) 



(59) 



/ JV 2 

- min Trjv/jy V |^ + ?7 + /T 1 In /a 



vi=l 



which we write as 
-Fo = min 



drw(r)/0Ar(r) + F N [p N 



(60) 



where the intrinsic Helmholtz free energy functional in 
the canonical ensemble is defined as 



N 



F N [p N }= min Tr N f N {J2g- + U + p- 1 \nfi f ], 

In^Pn J 

(61) 

which is formally equivalent to the definition (|2"2"|) of J 7 ^ in 
the grand ensemble upon identifying the different traces 
and different types of many-body distributions. It is 
clear that the density distribution pn.o that minimizes 
the RHS of (|60[) is the true equilibrium distribution in 
the canonical ensemble 



P/v,o(r) = Tr N fN,of>( r )- 



(62) 



and that Fq = Fjv [A/v,o] ■ The variational principle (foT))) 
implies that 



5F N [p 



Sp N (r) 



+ v(r) = 0, 



(63) 



where the derivative is taken under the constraint 
/ dr PN (r) = N. 

VII. DISCUSSION AND CONCLUSION 

The formulation of DFT rests on the existence and 
uniqueness of the intrinsic free energy as a functional of 
the one-body density for a given classical system. We 
have described two methods for defining this quantity, 
via Eq. (TT5l) based on the Mermin-Evans argument [3, 
7], and via (|2"2"|) based on Levy's constrained search [Tll |. 
Following the derivations presented in Sees. IIIII and IIV1 
it is clear that these methods are different in procedure 
and underlying principles. 

In the Mermin-Evans sequence of arguments it is for- 
mally proved that the equilibrium many-body probabil- 
ity distribution, /o is a functional of the equilibrium one- 
body density, po. The existence of this functional rests on 
a sequence of functional dependencies. For given inter- 
atomic potential U and given one-body density p, there 
is a unique external potential v, that generates this p. 
When input into the form of the many-body distribution 
in the grand ensemble (JIJ, this uniquely determines /o 
as used on the RHS of the definition (TTST) of the intrinsic 
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free energy functional J-[p]. This chain of dependency 
is implicit in order to properly define the free energy 
functional via (fT5|) . Note that the naive view that the 
equilibrium probability distribution, /o, is a function of 
the external potential, v, such that (|T5|) should also de- 
pend on v, gives the impression that functional (|15[) is not 
independent of the external potential energy. Certainly 
this is not the case-as one may recall the argument above 

On the other hand, Levy's method does not rely on 
the above rather subtle argument. An appealing feature 
of the constrained search method is the definition ([212]) of 
J-l[p]. Here the intrinsic free energy functional is explic- 
itly independent of the external potential, which is not 
as easily observed from T[p] of Eq. (fl~5j) . Kohn [l3| and 
Levy |18| describe the constrained search method as a 
two step minimization procedure, and we have laid out 
analogous reasoning in Sec. [Vj 

The underlying principle of the Mermin- Evans method 
of defining the intrinsic free energy functional is v- 
representability of the trial density, whereas Levy's 
functional is based on the weaker condition of /- 
representability. However, one may restrict the con- 
strained search to the class of one-body densities that 
is w-representable. In this case Levy's functional (|22|) 
becomes equal to the Mermin- Evans functional (JX5J . 
Hence, the constrained search method reduces to find- 
ing the equilibrium one-body density which correspond 
to an (equilibrium) external potential, v, that minimizes 



functional ([24]) over all one-body densities, p, each associ- 
ated to a specific v. Furthermore, applying the Legendre 
transform upon functional (|24[) and minimizing a set of 
external potentials which yields a fixed one-body density, 
gives Week's free energy functional [19]. 

In practice, minimizing Levy's version of the free en- 
ergy functional (|2"2"j) will certainly not be easier than 
solving the many-body problem itself. Hence, whether 
the definition (|22j) helps to construct approximations for 
grand-canonical free energy functionals remains an open 
question. For the case of the canonical ensemble we point 
the reader to the very significant body of work that has 
been carried out to formulate a computational scheme 
that permits to capture the effects that arise due to the 
constraint of fixed number of particles (2lT[26| . While the 
generalization to equilibrium mixtures is straightforward, 
we expect the application of Levy's method to DFT for 
quenched-annealed mixtures |28l - l3l1 ] to constitute an in- 
teresting topic for future work. 
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